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Abstract 

Recently a partially solvable many-body problem with nearest and next-nearest 
neighbour interactions is proposed JjJ]. We show that by adding a suitably chosen 
momentum dependent nearest neighbour interaction, such a model can be converted 
into an integrable system with Lax operator formulation and related conserved 
quantities. We also solve the eigenvalue problem for the model exactly and as a 
byproduct obtain some identities involving associated Laguerre polynomials. 

o ' 

It is well known that random matrix models have deep connection with some ex- 



actly solvable many-body systems having long-range interactions, namely Calogero 



a 

and Sutherland model 0. However recently a new type of random matrix model, known 

£h i 

as short-range Dyson model, has been introduced, where only the nearest neighbour en- 
ergy levels interact with each other |4], [|]. For finding out the corresponding many-body 
system, Calogero and Sutherland like models with nearest and next-nearest neighbour 
interactions have been proposed recently [|IJ. Though many interesting results, e.g. con- 
struction of static correlation functions, partial solution of the energy eigenvalue problem, 
existence of off-diagonal long-range order etc. have been obtained, the integrability and 
complete solvability of these models could not be established. The Hamiltonian of such 
Calogero like model as considered in |], |6| is given by 

1 N d 2 to 2 N ^ N ~ 1 1 N ~ 1 1 

2j =1 dx 2 2 - =1 3 - =l (xj — Xj + \) 2 - =2 — Xj)(xj — Xj + i) 

(1) 

Our aim here is to propose an integrable and completely solvable variant of the above 
model by suitably redefining its coupling constants and introducing an additional mo- 
mentum dependent nearest neighbour interaction. The Hamiltonian of this model may 
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be given as 



1 N pa , ,2 N N-l 



N-i y N d 

+ 92^2 / - _ r .w r ._ r , \ + ^ E/i flZ7> ( 2 ) 

where the last term represents a momentum dependent interaction with 

fj = (l-S v )—^— + (l-S NJ )—^— (3) 

and the coupling constants are parametrised as 

gi = v{% + u), g 2 = -v 2 , g 3 = -iv. (4) 

It is easy to see that for real parameter v the Hamiltonian (HD becomes hermitian. It 
has been shown recently JFj that the original Calogero model with long-range momentum 
dependent interaction retains its integrability and complete solvability. Similarly for the 
present model (|^), involving only nearest neighbour momentum dependent interaction, 
one can establish the integrability and the existence of infinite number of conserved quan- 
tities. For this purpose one needs the associated Lax operator, which may be expressed 
as a JV x JV matrix with elements 

d 

Ljk = ^~ i dx~ + U ^ 6jk ' ^ 
One may check that this operator and Hamiltonian (Q) satisfy the Lax equation 

[H, L ± ] = -LujL^, where L = Ljf. ± iuxjSjk- (6) 

Note that though the complementary Lax operator M is absent in this Lax equation, it still 
leads to a set of conserved quantities given by I m = J2j,k ((^ + ^~) m )jk > = 1, 2, 3, . . . , 
ensuring the integrability of the system. Evidently Hamiltonian (^) is the first element of 
this set. 

For solving the energy eigenvalue problem of model ([!]) and its other variants, the 
authors of |6| have adopted two different approaches. The first one, as described in 
[p[-|Jll|, constructs the eigenf unctions by mapping the model to a system of free oscil- 



lators through a suitable similarity transformation. In the second approach, pioneered 
by Calogero @], the solution for eigenfunctions is sought through an ansatz in the fac- 
torised form: Pfc(x)$(r 2 ), where Pk{x) is a homogeneous polynomial in all coordinates 
and $(r 2 ) is a function of r 2 = J2i=i %\- Through both these methods 
some particular but exact eigenvalue solutions were obtained in |], ||. 
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Due to the integrability property of our model (0), we hope to get the complete solution 
of its eigenvalue problem by adopting similar methods. Following the first approach, we 
find that there exists an unitary transformation which reduces the Hamiltonian @ to a 
system of free oscillators: 

'n-i \ iu {N-i \- iV 

llixj-Xj+i)] H f [] ( Xj -x j+1 ) J = H free , (7) 

with 

Hfree = ^E^ + yE^ (8) 
z j=1 UXj z j=1 

As is well known, the eigenfunctions for the free oscillator model may be given by 

N 

4>nx,n2,-,n N = ]J H n 3 ( V^Xj) , (9) 

where r 2 = Ylj=i x ] an d H n (x) is the Hermite polynomial of degree n. Therefore one 
can immediately write the eigenfunctions for our model (|2|) as 

( n - 1 V" ^ N 

\i=i / i=i 

The corresponding eigenvalues are evidently same as that of the free oscillators: 

(n n \ 

E ni>na ,-, nN = W ly + H™ij • (11) 

The ground state wave function is clearly given by ip gr = (j\f=i(xj — ^j+i)) e - ^, 
with energy eigenvalue E gr = Thus we see that model @, which is obtained from (P 
by adding a momentum dependent term and redefining the coupling constants, becomes 
exactly solvable giving very simple spectrum. In a similar way, by adding suitably chosen 
momentum dependent interactions, the other variants of model ([!]) considered in |T], [| 
can also be transformed to integrable systems and solved exactly. 

We should note here that the unitary transformation (|7|), which brings our Hamilto- 
nian to the free oscillator model, does not contain any singular terms. Consequently the 
eigenfunctions fllCf ) form a complete set and the spectrum coincides with that of distin- 
guishable oscillators. In contrast, the similarity transformations which map the standard 
Calogero model as well as ([!]) into the free oscillator model contain singular terms. In the 
case of Calogero model such singularities can be avoided by symmetrising the eigenfunc- 
tions of free oscillators in all coordinates ||. As a result the eigenfunctions of Calogero 
model represent a complete set with the excitation spectrum becoming equivalent to that 
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of bosonic oscillators. However in the case of (|1|) such singularities can not be removed 
in general even by taking symmetrised eigenfunctions and that makes the completeness 
of eigenfunctions difficult to establish. 

It may be noted that both ([I]) and our model (||) reduce to the free oscillator model 
at v — » 0. Therefore it may be expected that a link should exist between the solutions 
of these two models at v — > limit. To examine this possibility we follow the second 
approach, as mentioned above, for solving the model (@) and find that the corresponding 
eigenfunctions can be obtained in the factorised form 



,2 



'N-l , x 

-+k-l, 2 \ 



$n,fc = e-V I Y[(x j -x j+1 )j L r y (ur 2 )P k (x). (12) 
Here Pk{x) is a homogeneous polynomial of degree k satisfying the equation 

N Q2 

-K~2 P k(xi,X 2 ,- ■ ■ ,Xn) = 0, (13) 
3=1 ° X 3 

and Ln +k 1 (ur 2 ) is the associated Laguerre polynomial. The corresponding eigenvalue 
would be E n ^ = u + 2n + k). A class of trivial solutions of (fT3|) with k < N is 
clearly given by 

where ais are all different integers ranging from 1 to N. For finding nontrivial solutions 
of fll3D we express Pk{x) in a general form 

k\ 

Pk(x) = Yl i ! • • • r ! ^ ri ' T2 "-' rN x i lx 2 2 ' ' ■x'n i 

fi+^H hrjv=fc 1' 2 ' 

where f£ r2 ... rjv are yet unknown coefficients and the summation variables rj are nonneg- 
ative integers. The equation (|T3| ) however puts constraints on these coefficients as 

JV 

fr 1 ,r 2 ,--,r j - 1 , rj + 2,r j+1 ,r j+ 2, — ,r N ~ ^ ' 

(16) 

3=1 

where {rj} represents any partition of k— 2 satisfying J2f=i r j — k — 2. Such constraints are 
evidently different for different choices of {rj}. However since the number of constraints 
is less than the number of unknown coefficients, they can not be fixed uniquely unless one 
imposes some extra relationship among them. Therefore the corresponding polynomial 
(|T5D may not be symmetric in all its coordinates in general. For example for k = 2 one 
gets the solution 

P2(x) = aix\ + 02^2 + ■ • • + aNX 2 N , with dj = 0, (17) 

3=1 
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which can not be presented in a symmetric form. On the other hand for higher k one may 
obtain symmetric polynomials as solutions of (|16|). For k = 3 such a solution is explicitly 
given by 

N 

P 3 (x) = a$>? + &5>£c fcj (18) 
where 3a + (iV — 1)6 = 0. For k = 4 we get however two independent solutions of the form 

N 

P 4 (x) = a$>J + (19) 
j=i j^k 

with the constraint 3a + (N — l)b = 0, and 

P 4 (x) = c Y x f x k + d Y x)x k xi , (20) 

with 3c + (JV-2)d = 0. 

We may compare now the above solutions at v — > with those of ([!]), as obtained in 
[|IJ, at the same limit. For the case k = 2 while we get two solutions with P 2 (x) taking 
the trivial form ( |i~4"D as well as the nontrivial one flTFD, [Q| gives only one solution with 
trivial form of P2(x). For k = 3 we get again two independent solutions corresponding 
to (0) and (0), while in |I[ one finds only a particular combination of them. Similarly 
for k = 4 we get three independent solutions in contrast to a single solution of This 
picture prevails also for higher k. 

Finally we would like to comment on some useful identities, which may be derived by 
comparing the eigenf unctions (110 ) and (TT) obtained through two different approaches. 
Expressing each of the solutions (|12|) as some linear combinations of the complete set 
(|ToD , one would generate such an identity involving associated Laguerre and Hermite 
polynomials. By using this procedure along with the familiar conversion formulas like 
H 2n (x) = (-l) n 2 2n n\LnHx 2 ), H 2n+l (x) = (-l) n 2 2n+l n\xLl(x 2 ), we find interestingly 
that (O) with the trivial solution (EH) for Pk{x) leads to a well known summation formula 



12jl for associated Laguerre polynomials: 

V,, , k 1 N _1 

L* +k \tor 2 ) = Yl \{LUujx 2 r ) J] LnHux 2 r ). (21) 

n 1 +n 2 A hn N =n r=l r=k+l 

On the other hand, (0) with nontrivial solutions for Pk{x) leads to apparently new 
identities involving associated Laguerre polynomials. For example, the solution flT7| ) for 
N = 2 giving P2(x) = x\ — x\ results an identity of the form 

1 71+1 _i _i 

L 2 n {ux\ + ux 2 2 ) = — 2 E ( 2s " n ~ l)Lj s+l {ux\)L s *{ux 2 2 ) . (22) 

u\x x — x 2 ) s=0 

We hope to find similar identities for higher values of k and iV and report elsewhere [13|. 
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